A THIN FILM APPROXIMATION OF THE MUSKAT PROBLEM WITH 
GRAVITY AND CAPILLARY FORCES 



PHILIPPE LAURENQOT AND BOGDAN-VASILE MATIOC 

Abstract. Existence of nonnegative weak solutions is shown for a thin film approximation of 
the Muskat problem with gravity and capillary forces taken into account. The model describes the 
space-time evolution of the heights of the two fluid layers and is a fully coupled system of two fourth 
order degenerate parabolic equations. The existence proof relies on the fact that this system can 
be viewed as a gradient flow for the 2— Wasserstein distance in the space of probability measures 
with finite second moment. 



1. Introduction and main result 

The Muskat problem is a free boundary problem describing the motion of two immiscible fluids 
with different densities and viscosities in a porous medium (such as intrusion of water into oil). It 
gives the space and time evolution of the heights / > and g > of the two fluid layers, with the 
first layer, of height /, located on a impermeable horizontal bottom and the second one, of height 
g, on top of it, as well as that of the pressure fields inside the fluids. As it involves four unknowns 
and two free boundaries, one separating the lower and the upper fluid and one separating the upper 
fluid and the air, it is a complex problem. Recently, using a lubrication approximation, see, e.g., 
[14, Chapter 5.B], a thin film approximation to the Muskat problem has been derived in [10] which 
retains only the heights / and g of the two fluid layers as unknowns and reads 

f d t f =div [/ (-AVAf - BVAg + aVf + bVg)} , 

<^ (t, x) e (0, oo) x Mr, (1.1a) 

I d t g =div [g (-VA/ - VAg + cV/ + cVg)) , 

together with initial conditions 

(/,3)(0) = (/o,3o), x£R 2 . (Lib) 

The parameters (A, B,a, b, c) involved in (1.1a) depend on the densities, viscosities, and surface 
tensions of the fluids and are assumed to satisfy 

(a, b, c) E [0, oo) 3 , cB = b, and A > B > 0. (1.2) 

Let us recall that the second order terms in (1.1a) account for gravity forces while the fourth order 
terms result from capillary forces in the original Muskat problem. 

The problem (1.1a) is a fourth order degenerate parabolic system with a full diffusion matrix. Its 
parabolicity has been exploited in [11] to show the local existence and uniqueness of positive strong 
solutions to (1.1a) in a bounded interval (0,L) with no slip boundary conditions. Global existence 
for initial data close to a positive flat steady state is also proved when ir 2 (A — B) jl? + (a — b) > as 



Date: March 8, 2013. 

2010 Mathematics Subject Classification. 35K65, 35K41, 47J30, 35Q35. 

Key words and phrases, thin film, degenerate parabolic system, gradient flow, Wasserstein distance. 

1 



2 



PH. LAURENQOT AND B.-V. MATIOC 



well as the stability of this steady state. The existence and stability of flat and non-flat stationary 
solutions are also discussed according to the values of the parameters. Nonnegative global solutions 
have also been constructed either when only gravity is taken into account (A = B = 0) and (1.1) 
is considered in a bounded interval with homogeneous Neumann boundary conditions [9] or in M 
[13], or when gravity forces are discarded (a = b = c = 0) and (1.1) is considered in a bounded 
interval with no slip boundary conditions [16]. An important tool in the above mentioned works is 
the availability of two Liapunov functionals for (1.1), one being the functional E defined by 

£(f,9)-=l [ [(A-B)\Vf\ 2 + B\V(f + g)\ 2 + (a-b)f 2 + b(f + g) 2 ] dx (1.3) 

for (/, g) € i7 1 (M 2 ;IR 2 ) which decreases along the trajectories of (1.1) according to 

±-8{f,g) = - [ \f\VA(Af + Bg)-V(af + bg)\ 2 + Bg\VA(f + g)-bV(f + g)\ 2 
at j R 2 l 



dx . 



It turns out that there is an underlying structure in (1.1) which allows us to view it as a gradient 
flow for the energy £ defined in (1.3) with respect to the 2— Wasserstein distance in V 2 (R 2 ) x V 2 (M. 2 ). 
Recall that V2CR 2 ) is the set of Borel probability measures on IR 2 with finite second moment and 
that, given two Borel probability measures \i and v in V 2 ^ 2 \ the 2— Wasserstein distance W<i{p,, u) 
on V 2 (R 2 ) is defined by 



W£((J,,v):= inf / \x - y\ 2 d7r(x, y) , 
7ren(/i,i/) Jr4 



where II(/x, v) is the set of all probability measures tt £ V(R ) which have marginals /i and u, that 
is, ir[U x R 2 ] = fj,[U] and vr[lR 2 X U] = v[U] for all measurable subsets U of R 2 . 

This gradient flow structure was actually uncovered in [13] in the simplified case where the 
capillary forces are neglected (A = B = 0) and shown to provide a convenient setting to construct 
weak solutions to (1.1). This is thus the approach we shall use in this paper to prove the existence 
of nonnegative weak solutions to (1.1), showing additionally the convergence of the variational 
approximation. It is worth noticing at this point that, if g = 0, the system (1.1) reduces to the 
single equation 

d t f = div [/ (-AVAf + aV/)] , 

which also has a gradient flow structure with respect to the 2— Wasserstein distance as already 
shown in [18, 19] for A = and [17] for a = 0. Let us recall that, since the pioneering works 
[12, 18, 19], several parabolic equations have been interpreted as gradient flows for Wasserstein 
distances, including the Fokker-Planck equation [12], the porous medium equation [19], second order 
nonlocal and/or degenerate parabolic equations [1, 3, 4], some kinetic equations [6, 8], and fourth 
order degenerate parabolic equations [15, 17]. Besides (1.1) there does not seem to be many systems 
of parabolic partial differential equations endowed with a similar structure with the exception of 
the parabolic-parabolic chemotaxis Keller-Segel system which has a mixed Wasserstein-L2 gradient 
flow structure [5, 7]. 

Before stating the main result of this paper, we introduce some notations: let 1C be the convex 
subset of V 2 {R 2 ) defined by 

K := {h € Li(R 2 , (1 + \x\ 2 )dx) n H l {R 2 ) : h > a.e. and \\h\\i = l} , 
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and set /C 2 := K, x /C. We next recall that the functional H defined by 

H(h) := / h\n(h)dx (1.4) 

is well-defined for h G /C, see Lemma A. 3. 

The main result of this paper is the following: 

Theorem 1.1. Assume £/ia£ (1-2) is satisfied and let (fo,go) G K-2- Given r G (0,1), we define 
(f®,g®) := (fo>9o) an d consider for each n G N £/ie minimization problem 

inf J?(u,i;), (1.5) 

u>/iere 

J7(«,«) := i- [Wf(u, + £ Wf («,<£)] + £(«,«), («,«) G/C 2 . 

For eac/i n G N, £/ie minimization problem (1.5) /ios a solution (/" , <?™ +1 ), which is unique if 
a > b. Defining the interpolation functions (/ T , <7t) (/t> 5t)(£) := (fridr) f or t ^ [ nr > ( n + l) r ) 
and n G N, £/iere e:cis£ a sequence \ and functions (/,<?) ■ [0, oo) — >• /C 2 suca £/ia£ 

(f Tk ,gr k )(t)^(f,g)(t) m L 2 (M. 2 ;R 2 ) (1.6) 

/or a// £ > 0. Moreover, 

(*) (/,5) GL oo (0,£;i/ 1 (K 2 ;^ 2 ))nL 2 (0,£;^ 2 (]R 2 ;]R 2 )), 
(ii) (/,<?) G C([0,oo),L 2 (M 2 ;lR 2 )) and (/, 5 )(0) = (/ , 5o ), 
and £/ie pair (7,5) is a u>eaA; solution of (1.1) in £/ie sense that 



I (/(*) - /o)£ dx+ f [ [A(Af + B 5 ) div(/V£) + /V(a/ + 6a) • V£] dx ds = 

Jr 2 jo jr 2 



1.7) 



/ / [A(/ + 5 ) div( o V£) + caV(/ + a)-V£]dxds = 

JO JR 2 

/or all t > and £ G C X) (M 2 ). Finally, (f,g) satisfies the following estimates: 

H(f(t)) + £?#(<?(£)) + T D H (f(8),g( a )) ds < H(f ) + 5#( go ) , (1.8) 

Jo 

£(f(t),g(t)) + \ f (IK(s)H! + B\\w g {s)\\l) ds < £(fo,g ) (1.9) 
/or a// £ > 0, where 

D H (f,g) := (4 - £)||A/|| 2 + B||A(/ + ff)||| + (a - 6)||V/|| 2 + 6|| V(/ + g)f 2 , 

and and w g are two vector fields in iv 2 ((0, 00) x R 2 ; M 2 ) defined as follows: introducing the vector 
fields (jf,j g ) defined by 

j f := -V(/A(A/ + B 5 )) + A(Af + Bg)Vf + fV(af + bg) 

in V'((0, 00) x M 2 ;IR 2 ), (1.10) 

i 9 := -V(aA(/ + a)) + A(/ + g)Vg + caV(/ + g) 
they actually both belong to L 2 (0, 00; L4/3(IR 2 ; M 2 )) and 

jf = \/~f wj and j g = ^/g w g a.e. in (0,oo) x f 2 . (1-H) 
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Notice that Theorem 1.1 provides not only the existence of a weak solution to (1.1) but also the 
convergence of subsequences of solutions to the variational scheme (1.5) to a solution to (1.1). The 
proof of Theorem 1.1 thus relies strongly on the study of the minimization problem (1.5) which is 
performed in Section 2. The availability of the second Liapunov functional H(f) + BH(g) (which 
is really a Liapunov functional only if a > b) plays an important role here as it allows us to show 
by an argument from [17] that the minimizers of (1.5) are actually in H 2 (M?). This additional 
regularity is actually at the heart of the identification of the Euler-Lagrange equation satisfied by 
the minimizers, see Section 3. The convergence of the scheme and the existence of weak solutions to 
(1.1) are established in the last section, some care being needed to handle the fourth order terms. 
Indeed, as indicated in (1.10) and (1.11), due to the degenerate character of the equations, we have 
to deal with expressions that correspond to the fourth order terms in (1.1a) and are not functions. 

Remark 1.2. (1) In contrast to the dissipation of the energy £, the dissipation Djj(f,g) of the 
functional H(f) + BH(g) need not be nonnegative and this occurs when a < b, see (1.8). It 
is yet unclear what information on the dynamics may be retrieved from (1.8) in that case. 

(2) For simplicity, we have restricted the analysis to initial data (/o,5o) satisfying ||/o||i = 
||<7o||l = 1- A similar result holds true for arbitrary nonnegative and integrable initial 
data (/o,5o) £ i^ 1 (K 2 ;M 2 ) with finite second moment and may be proved analogously to 
Theorem 1.1. Introducing (F,G) := (//||/o||i, ff/Hffolli)) this new unknown solves a system 
with the same structure as (1.1), but with different parameters (depending on ||/o||i and 

Ibolli)- 

(3) Theorem 1.1 is also valid for the one dimensional version of (1.1). 

Throughout the paper we use the following notation: x = (xi,X2) denotes a generic point of 
M 2 and the partial derivative with respect to X{ is denoted by <9j, i = 1,2. For a sufficiently 
smooth function h, V/i := (d±h,d2h) denotes its gradient and D 2 h := (9j5j/t) 1<i <2 its Hessian 
matrix. Finally, D£ denotes the gradient of a vector field £ = (£1,^2) £ C 1 (IR 2 ;R 2 ) and is given by 
D£ := (<9i£j) 1< ,- - <2 , the divergence of £ being the trace of D^, that is, div(£) = <9i£i + c?2^2 ■ 

2. The minimization problem 

In this section, we show that, given (/o,<7o) & ^2 an d T > 0, the minimization problem 

inf T T (u,v) (2.1) 
(u,v)e!C2 

with the functional T T defined by 

T T (u,v):=^ (W^uJ ) + BW^v,g ))+S(u,v), (u,«)e/C 3 , (2.2) 

has at least one solution and we study the regularity of the minimizers. Since a — b might be 
negative, which is the case when the more dense fluid lies on top of the less dense one, cf. [11], the 
first step is to prove that £ is bounded from below. 

Lemma 2.1. The energy functional £ satisfies 

£(/,$) >-Ci + ! HV(/ + 5)ll2 + ^-^ HV/Hi forall(f,g)elC 2 , (2.3) 
where C\ is a positive constant depending only A, B, a, b, and c. 
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Proof. According to the Gagliardo-Nirenberg-Sobolev inequality there is a positive constant C 2 > 
such that 

\\h\\ 2 < C 2 \\h\\\ /2 \\Vh\\y 2 for all h G ^(R 2 ). (2.4) 
By the definition of /C, this inequality yields 

IHI2 < C l ll V/i ll2 for all h€K. (2.5) 
Owing to (2.5) and Young's inequality, we find, for (/, g) G IC 2 , 

2£(f, g) = (A - B) \\Vfg + B ||V(/ + g)\\l + (a - b) \\f§ + b\\f + gg 
>B ||V(/ + g)\\l + {A-B) || Vf\\l - C\(b - a)+ ||V/|| 2 

>B ||V(/ + 5 )||I + ^^ ||V/|| 2 

, A-B / C|(6-a) + V Cf(6-o) 



2 V A-B 7 2(A-B) ' 

whence (2.3). □ 

We now show the existence of a minimizer to (2.1). 

Lemma 2.2. Given (/o,<7o) G ^2 anc ^ t > 0, f/iere exists a minimizer (f,g) G /C2 0/ (2.1) which is 
unique if a> b. Moreover, f and g both belong to H 2 (R 2 ) and 

(A-B) \\Afg + B \\A(f + g)g + (a - b) ||V/|| 2 + b\\V(f + g)g 
< i [H(f ) - H(f) + B (H(g ) - H(g))} , (2 ' 6) 
the functional H being defined in (1.4). 

Proof. Pick a minimizing sequence (u k ,Vk)k>i of J> in /C2- Invoking (2.3) and (2.5) we see that 

\\ukl\H 1 + IKIIifi < C, k>l, (2.7) 
W 2 (u k J ) + W 2 (v k ,g ) < C, k>l, (2.8) 
and the estimate (2.8) implies that 

(u k + v k ){x){l + \x\ 2 ) dx < C , k>l, (2.9) 

by classical properties of the Wasserstein distance W 2 , see, e.g., [12, Eq. (17)]. The compactness of 
the embedding of H 1 (R 2 )nLi(R 2 ; \x\ 2 dx) in L 2 (R 2 ) (which is recalled in Lemma A.l below) ensures 
that there are non- negative functions / and g in H 1 (R 2 ) n Li(R 2 , (1 + |x| 2 ) dx) and a subsequence 
of (u k ,v k ) k >i (not relabeled) having the property that 

(u k ,v k ) -> (f,g) in L 2 (R 2 ;R 2 ) and a.e. in R 2 , 

(2.10) 

(u k ,v k )^(f,g) mH l (R 2 ;R 2 ). 

Furthermore, it readily follows from (2.7), (2.9), and the Dunford-Pettis theorem that (u k ) k >\ and 
{vk)k>i are weakly sequentially compact in L\(M. 2 ). This property, together with (2.10) and the 
Vitali theorem imply the strong convergence of (u k )k>i anci {v k ) k >i m Li(R 2 ), namely, 

(u k ,v k )^(f,g) inL!(IR 2 ;M 2 ). (2.11) 
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Since (uk,Vk)k>i belongs to IC2 for each k > 1, we conclude from (2.9). (2.10), and (2.11) that 
(/iff) ^ ^-2- That (f,g) is a minimizer of JF T in /C2 follows from (2.10), (2.11), the weak lower 
semicontinuity of S, and the fact that the 2-Wasserstein metric W2 is lower semicontinuous with 
respect to the narrow convergence of probability measures in each of its arguments, the latter 
convergence being guaranteed by (2.11). This completes the proof of the existence of a minimizer 
to T T in K,2- 

Next, when a > b, the uniqueness of the minimizer follows from the convexity of JC 2 and W 2 and 
the strict convexity of £. 

We finish off the proof by showing that any minimizer (/, g) of T T in K 2 actually belongs to 
if 2 (M 2 ;M 2 ), the proof relying on a technique developed in [17]. More precisely, denoting the heat 
semigroup by (G s ) s >q which is defined by 



(G s h)(x) := J exp A f ) %) dy , (a, x) G [0, 00) x 



d2 



for h G Li(R 2 ), classical properties of the heat semigroup ensure that (G s f, G s g) G 1C 2 for all s > 
since (/, g) G K-2- Consequently, J- T {f, g) < J- T {G s f ,G s g) and we deduce that 

£(f,g) - £(G s f,G s g) < i- [(W|(G S /, / ) - Wf (/,/ )) + B (W^G^So) - ^ 2 (<?, 5o ))] (2.12) 

for all s > 0. On the one hand, an explicit computation gives for s > 
d 



ds 



£(G s f, G s g) = [ [(A- B) VG s f d s (VG s f) + B VG s (f + g) d s (VG s (f + g)) 
+(0 - 6) G s / a s G s / + 6 G s (/ + 5 ) d s G s (f + 5 )] d.x 



(A -5) ||AG s /|| 2 -i? ||AG S (/+- 



12 

-(a -6) ||VG S /|||-6||VG S (/ + 5 )|||, 
which yields, after integration with respect to time, 

- f [(A - B) \\AG a f\\ 2 2 + B \\AG a (f + 9)\\i + (a - b) ||VG CT /|| 2 + b \\VG a (f + g)\\ 2 2 ] da 
s Jo 

= £(f,g)-£(GJ,G s g) 
s 

Since a 1— >• || ACr CT /i||2 and c 1— > ||VG CT /i||2 are non-increasing functions for all h G Li(R 2 ), we end up 
with 

{A - B) ||AG S /|| 2 + B ||AG S (/ + 5 )|| 2 + (a - 6)+ ||VG S /|| 2 + b \\VG s (f + 5 )|| 2 

£(f,g)-£(G s f,G s g) „ . ( 2 -!3) 
< w.y; 1 + (6 - a )+ V/ I 

s 

for all s > 0. On the other hand, since the heat equation is the gradient flow of the entropy functional 
H defined by (1.4) for the 2-Wasserstein distance W2, see, e.g., [2, 12, 19, 21], it follows from [2, 
Theorem 11.1.4] that, for all (h,h) G K 2 , 

1 Wi(G s h, h) + H(G s h) < H(h) for a.e. s > 0. (2.14) 

2 ds 
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With the choices (h,h) = (/, /o) and (h,h) = (g,go) in (2-14), we obtain 

\js [ W i( G sfJo) + BW^(G s g,g )] < H(f ) - H(G s f) + B (H(g ) - H(G s g)) 

for a.e. s > 0. Integrating the above inequality with respect to time and using the time monotonicity 
of s i — y H(G s f) and s *— > H(G s g) lead us to 



1 

2s 



[Wi(G s f, /o) - Wi(f, f ) + B (W 2 (G s g,g ) - W 2 (g,g ))} 



< [H(f ) - H(G s f) + B (H(g ) - H(G s g))] 
Combining (2.12), (2.13), and (2.15), we find 



(2.15) 



(A-B) \\/\G s f\\l + B \\AG s {f 



+ (a-6)+ ||VG,/||| + 6 ||VG S (/ + , 



< - [Htfo) - H(G s f) + B (H(g ) - H(G a g))] + (b - a)+ ||V/||| 



(2.16) 



for s > 0. Since (/, g) G /C2, classical properties of the heat semigroup and the functional H entail 
that (H(G s f), H(G s g)) converges towards (H(f), H(g)) as s — > 0. This convergence, (2.16), and 
the positivity of A — B and B readily imply that (AG s /) s >o and (AG s (f + g)) s >o are bounded in 
L2(M 2 ). Since they converge towards Af and A(/ + g), respectively, in the sense of distributions 
as s — > 0, we deduce that both A/ and A(/ + g) belong to L2(M?) and (AG s f, AG s (f + g)) s >o 



converges weakly to (A/, A(/ + g)) in 



p2. up 2^ 



as s — > 0. As a consequence, recalling that 



(f)9) G /C2, we conclude that / and f + g belong to ff 2 (]R 2 ), and so does g. It remains to pass to 
the limit s ->■ in (2.16) to obtain (2.6). □ 

3. The Euler-Lagrange equations 

We now identify the Euler-Lagrange equation satisfied by the minimizers of the functional J- T , 
defined in (2.2), in A^- 

Lemma 3.1. Consider (/o,<?o) £ K-2 an d r > 0. If (f,g) is a minimizer of T T in JC2, it satisfies 
[ {f-fo)£dx + T f [A(Af + Bg) div(/V£) + / V(of + bg) ■ V£] fix 

i^ 2 eii 



(3.1) 



< 



anc? 



/ (g - 5o )e + r / [A(/ + g) div( 5 V^) + c<? V(/ + </) • V£] dx 

/ ll^ll 



(3.2) 



^2(5,50) 



/or £ G C^°(M 2 ), where D 2 ^ denotes the Hessian matrix of £. 



Proof. By Brenier's theorem [21, Theorem 2.12], there are two measurable functions T 
and 5 : E 2 -»• M 2 such that (/,#) = (T#/ , 5#ff ) and 



^ 2 2 (/,/o) 



k - r (x)| 2 fo(x)dx and W 2 2 (0,0o) 



|x — <S(x)| 2 go(x)dx . 



(3.3) 
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We pick now two test functions r\ = (771,772) and £ = (£1,^2) in Cg°(R 2 ;R 2 ) and define 
(T e , S e ) := (id +s:£, id +£77) , 

(f £ ,g £ ) := ((T e oT)#/o,(S £ oS)#s ) = (T £ #f,S £ #g) , (3 ' 4) 

for each e £ [0, 1], where id is the identity function on R 2 . Clearly, there is £0 £ (0, 1) small enough 
(depending on both £ and 77) such that, for e G [0, £0], T e and S £ are C°°— diffeomorphisms from R 2 
onto R 2 with positive Jacobi determinants 



4 := det{DT £ ) = 1 + £ div(f ) + £ 2 det(Df ) , 
:= det(L>S £ ) = 1 + £ div(r?) + £ 2 det(£>? ? ) . 

By (3.4), we have the identities 



(3.5) 



fe= f e ° Te \ and g £ = f t , eG(0,£ ], (3.6) 
J| o T £ Je° 



from which we deduce that / e and g e both belong to (R ) and satisfy ||/ e ||i = ||/||i = ||ffe||i = 
||g||i = 1. Additionally, we compute 

ll/ £ ||l=/ ^-dx and \\g £ \\ 2 2 = f ^r\dx. (3.7) 

For further use, we now study the behaviour of (f £ ) £ and (g £ ) £ as e — > 0. To begin with, we 
notice that, given a test function # E Cg°(R 2 ), it follows from (3.5) and (3.6) that 



lim / f £ ti dx = lim / / (■& o %) dx = I f 



fidx (3. 



and 



r f — f r ■& oT — i? 
lim / - — i- tidx = ]im f dx = I fV$-£dx, (3.9) 

e->0/ ltt 2 £ e-*-0/ TO 2 £ 



while (3.5), (3.7), and the Lebesgue dominated convergence theorem entail that 

lim ||/ £ ||! Hl/H 2 . (3.10) 

Thanks to (3.8) and (3.10) on the one hand and to (3.9) and Lemma A. 2 on the other hand, we 
conclude that 

fe^f in L 2 (R 2 ) and ^ £ ~ - -div(/£) in L 2 (R 2 ) . (3.11) 



Similarly, we get 



# e -> g in L 2 (R 2 ) and — — - ->> -div( OT ) in L 2 (R 2 ) . (3.12) 



We next turn to the convergence properties of (V/ e ) e and (X7g £ ) £ . Differentiating (3.6) and using 
(3.5), we find 

Vf £ o T £ = -L- Vf + e V £ (f, - £ 2 -L- R £ (0 , (3.13) 
(J?) 2 (Je) 3 
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where V e (f,0 := (y 1)£ (/, £), V 2 , £ {f, £)), 2^(0 ■= (iMO>-*MO)> 

V/-(<9 2 &,-di£ 2 ) /^div(0 



(J 



^2 



(J 



^3 



V/- (-^1,5^1) /5 2 div(0 



(J 



I) 2 



(J 



I) 3 



and 



i?i, £ (£) := didet(L>£) + #idiv(£) d 2 & - d 2 div(£) 9i£ 2 + e (^det^) d 2 £ 2 - ^det(D^) d^ 2 ) 
R2AO ■■= &det(£>£) + £> 2 div(£) d^i - Sidiv(0 + e @,det(I>0 ftft - ftdet(Df) flbfi) 
Let us now draw several consequences of (3.13): first, we have 



liv/«||l = 

— > 1 in L 



|V/ £ 



oT e | 2 j| dx, 



and, since / G ^(IR 2 ) and J £ 
previous identity that 

lim||V/ £ ||| = ||V/||i. 

In particular, (V/ £ ) £ is bounded in L 2 (R 2 ;IR 2 ) and, since / £ — >■ / in L 2 
that (V/ £ ) £ converges weakly towards V/ in L 2 (M 2 ; M 2 ) . This convergence and (3.14) then guarantee 
that 



by (3.5), it readily follows from (3.13) and the 

(3.14) 

2 ) by (3.11), we conclude 



V/ e ^V/ in L 2 



(3.15) 



Next, owing to (3.13), we have 



1 



V(/ e - /) 



1 - Jl oT; 



+ 



+ V e (f,0°T- 1 -e 



f°T~ l 



1 



V/oT; 



Jl oTf 1 

f ^(e)oT; 1 . 



v/ 



The properties of jf, T £ , /, and the definitions of V s (f,£) and i? £ (£) readily ensure that the first, 
third and fourth terms on the right-hand side of the above identity are bounded in L 2 (R 2 ;R 2 ) while 
the boundedness in L 2 (M 2 ;K 2 ) of the second one follows from Lemma A. 2 since / £ H 2 (M?) by 
Lemma 2.2. Consequently, (V(/ £ — f)/s) e is bounded in L 2 (R 2 ;]R 2 ) and, recalling that ((/ £ — /)/e) £ 
converges weakly to — div(/£) in L 2 (M 2 ;M 2 ) by (3.11), we conclude that 

V(/ £ -/) 



Similar computations give 

Vg £ ->■ Vg in L 2 { 



and 



Vdiv(/0 in L 2 
V(<? £ - 9) 



')■ 



X7div(gt]) in L 2 



(3.16) 



(3.17) 



After this preparation, we can start the proof of (3.1) and (3.2). Recalling that (f e ,g s ) £ ^2 for 
all e 6 (0,£o], the minimizing property of (/, g) entails that J : T {f,g) < J~ r {.fei9e)i that is, 

0< ^ [W 2 (f e ,f )-Wi(fJ ) + B(W 2 (g £ ,g )-W 2 (g,g ))] + £(f e , g £ ) - 8 (/, g). (3.18) 
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Since (T e o T)#/ = f e by (3.4), we infer from (3.3) that 

W2(f e , f ) < WHf, f ) - 2e f (id — T) • (£ o T) f dx + e 2 [ \£o T\ 2 f dx. 

JR JR 
A similar inequality being valid with (g e , go, g, S £ , S,rj) instead of (f e , Jq, /, T e , T, £), we conclude 
that 

lim sup (1 [W 2 (f £ , f ) - W 2 (f, f ) + B {W 2 (g £ ,g ) - W 2 (g, go))] 



< 

We next show that 



/ (id-T)-(£oT) fodx + B [ (id -S) -( v oS)g dx 
Jr. 2 Jr 2 



lim 1 [(A - B) ||V/ £ ||I + fl || V(/ e + fc )||jj — {A — B) ||V/||| - S ||V(/ + ff )|||] 



/ [(AAf + BAg) div(/0 + BA(/ + (?) div( OT )] dx. 
Jr 2 



(3.19) 



(3.20) 



lim f [(a -b)\\f E \\l + b \\f £ + 5e ||l - ( a - 6) (I/HI - 6 ||/ + g\\j] 

f \ f 2 g 2 1 

= - / a — div(£) + 6 — div(r?) + 6 / div(g7/) + b g div(/£) . 

Jr 2 L 2 2 

Indeed, we write 

(a - b) \\f £ \\ 2 + b \\f £ + g £ \\ 2 - (a - b) \\f\\ 2 - b \\f + s||| = a Jf + b If + 6 If, 

with 

^:=ll/ £ ||!-||/ll!, ^:=W1-Nll, Jf := 2 /" (/ £ g e - f g) (x) dx. 

Jr 2 

It readily follows from (3.11) that 

lim f = lim / dx = - [ f div(/0 dx = -\f f 2 div(£) dx . (3.21) 

Similarly, (3.12) guarantees that 

lim ^ = lim / dx = -- ! g 2 div(r?) dx. (3.22) 

We next write 7| as 

/ [(fe + f)(9e-9) + (ge+g)(fe-f)]dx, 

Jr 2 

and deduce from (3.11) and (3.12) that 

lim §■ = - f [f div( OT ) + g div(/£)] dx. (3.23) 

Combining (3.21), (3.22), and (3.23) gives the claim (3.20). 
Finally, we show that 



(3.24) 



To this end, we write 

(A-B) ||V/ e ||! + fl \\V(f e +g e )\\1-(A-B) ||V/||| - B ||V(/ + g)f 2 = A L\ + B L\ + B L§, 



A THIN FILM APPROXIMATION OF THE MUSKAT PROBLEM 11 

with 

£iH|V/ £ ||i-||V/||| 5 £I:=||V<? £ |||-||V<?||L L e 3 :=2 [ (V f e ■ Vg £ - V / • Vg) dx. 

Jm. 2 

Thanks to (3.15) and (3.16), we have 



i im a =lim / v(/+/j.y(/ e -/) tiI 

e->o 2e e->o ha2 2 e 



f V/-Vdiv(/£)dx = / A/div(/£)dz, 
Jm. 2 Jm. 2 



(3.25) 



and similarly, by (3.17), 



lim ^ = lim I Yk±lA . Ykl^A dx= f Ag div(OT) dx . (3.26) 
e ^o2e e^oJ R2 2 e Jm 2 

Finally, 

L%= f Mfe + f) ■ Vfe - 5) + V(/ e - /) • V(g e + g)\ dx , 
Jm 2 

and we infer from (3.15), (3.16), and (3.17) that 

lim ^ = - f [V/ • Vdiv( OT ) + Vdiv(/£) • Vp] 

7k2 (3.27) 
[A/ div( OT ) + Ag div(/£)] dx. 

Combining (3.25), (3.26), and (3.27) gives the claim (3.24). 

We now divide (3.18) by e and take the limsup as e — > 0, using (3.19), (3.20), and (3.24). The 
resulting inequality being also valid for (—£,—7/), we obtain, by choosing successively £ = and 
77 = 0, that 

£, °T fodx = [ [A(Af + B<?) div(/£) + / V(a/ + 6 5 ) ■ £] dx , (3.28) 

r? o S g dx = [ [A(f + g) div( 5 £) + c / V(/ + 5 ) • £] dx . (3.29) 

Consider finally S G C^°(M 2 ) and take £ = V~ in (3.28). Since 

|S(x) - E(T(x)) - VS(T(x)) • (x - T(x))| < H^IU k~T(x)| 2 

for all x G ffi 2 by the mean value theorem, we find after multiplying the above relation by fo and 
thereafter integrating over K 2 and using (3.3) that 





R 2 






f (id 

R 2 


I) 



[~(x) - S(T(x)) - VH(T(x)) • (x - T(x))] / (x) dx 



Combining the above inequality with (3.28) gives (3.1). The inequality (3.2) next follows from 
(3.29) in a similar way. □ 
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In order to present the next result, we introduce first some notations. Given a nonnegative and 

,h 
5 



continuous function h and 5 > 0, we define the open sets V$ and V h by 



V h & := {xeR 2 : h(x) > 5} and V h := [j V h 5 . 

Lemma 3.2. Given (fo,go) G ^2 an d r > 0, any minimizer (f,g) of T T in IC2 is such that 
Af + Bg G Hf 0C (Vf ) and f + g G Hf 0C (V 9 ). Moreover, the functions jt, wj, j g , and w g defined by 

W f-\ a.e. R 2 \ 7>/ , ' ^ - V/ . (3-30) 

and 

w ._/ (-VA(/ + 5 ) + cV(/ + 5 )) a.e. zn , 

belong to L2(R 2 ;R 2 ) and satisfy 

[ [A(Af + Bp) div(/0 + / V(a/ + &</) -£]dx= I j f ■ £ dx, (3.32) 
Jk 2 Jm 2 

/ [A(/ + <?) divfoO + cs V(/ + (/)•£] dx= / j ff -e<te, (3.33) 

for all £ G C^°(]R 2 ; M 2 ). In addition, we have the following estimates: 

T\\wf\\ 2 <W 2 (fJo) and r|K|| 2 <^ 2 (<?,3o)- (3.34) 

Proof. Since iT 2 (]R 2 ) is embedded in C(R 2 ), Lemma 2.2 guarantees that (f,g) G C(R 2 ;IR 2 ) so that 
V f and pff are indeed open subsets of R 2 . Next, recalling (3.28), we use once more the embedding 
of H 2 (R 2 ) in C(R 2 ) as well as (3.3) to obtain that, for f G Cq°(]R 2 ; R 2 ), 



[ [A(Af + Bg) div(/0 + / V(o/ + bg) ■ £] dx 
Jm 2 

< \ (J I id ~ T \ 2 fo dx\ ' (J |e o T\ 2 f dx\ ' 

<^M([ |eiv^^<cW^ii/eiieii, 



We may thus extend the functional 

e — ► / [A(A/ + 5 5 ) div(/0 + / V(o/ + 65) • f] dz 
Jr 2 

to a continuous linear functional on L2(K 2 ;M 2 ). Consequently, there exists a unique function jt G 
L2(M 2 ;M 2 ) having the property that 

/ [A(Af + Bg) div(/0 + / V(of + &<?) • £] dx = / j/-^ for all £ G (^(M 2 ; R 2 ). (3.35) 
Jm 2 jr 2 

Since (/, g) G 7T 2 (]R 2 ;M 2 ) by Lemma 2.2, a density argument ensures that the relation (3.35) is 
actually true for all £ G i^QR 2 ; M 2 ). 



A THIN FILM APPROXIMATION OF THE MUSKAT PROBLEM 

Consider now 5 > and E G C^(V^;R 2 ). Clearly S/f € ^(R 2 ;^ 2 
with £ = E// that 
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/ A(Af + Bg) div(S) 



dx 



L 2 (P 



\V(af + bg)\\ 2 + \\E\ 



) and we infer from (3.35) 



(3.36) 



A duality argument then gives that A(Af + Bg) G H X (T>1) for all 5 > 0. Consequently, we get that 
Af + Bg £ Hf oc {T> f ) and together with (3.35) we deduce 



Jf 



-fVA(Af + Bg) + fV(af + bg) a.e. in . 



(3.37) 



We next prove (3.34), adapting an argument from [18, Proposition 2] and [13, Corollary 2.3]. Let 
X G Cq°(R 2 ) be a non-negative function with ||x||l = 1 an d set Xm( x ) := m 2 x( rn 'X) for i£R 2 and 
m > 1. Since H 2 (R 2 ) is embedded in C(]R 2 ), we have 



y m -=- + \\xm*f-f\\ 



1/2 



m 

Given G C^°(]R 2 ; M 2 ), the vector field •& / y/Y m + Xm * f belongs to C^(]R 2 ;M 2 ) too, and, by (3.3), 
(3.28), and (3.35) with the choice £ = "&j \JY m + Xm * 7i 



as m — )• oo . 

ioo/Tn)2. Tip 2> 



(3.38) 



if* 



V Y m + Xm* f 



dx 



< 



< 



W 2 (fJ ] 



+ Xm * / 



1/2 



W 2 (f,fo) 



f 



Ym + Xm * f 



1/2 



2- 



A duality argument then ensures that, for each m > 1, J/ / \/^m + Xm * / belongs to L2(M 2 ;M 2 ) 
with the estimate 



J/ 



y/Ym + Xm* f 



< 



W2UJ0) 



f 



Y m 4" Xm * f 



1/2 



Observing that 
< 



/ 



Y m + Xm * f 

we actually have the estimate 



f -Xm* f + Xm* f , 11/ ~ Xm * f\\oo -. , -. v 

^ T 1 i if J-m > 



Y m + Xm * f 

Jf 



Y 
1 it 



y/Y m + Xm * / 



<«(i t y,). 



(3.39) 



Several consequences can be drawn from (3.39): first, since Y m — > as m — > 00 by (3.38), 
the sequence (jf / \JY m + Xm * f)m is bounded in L2(R 2 ;M 2 ) and there are thus a subsequence 
of (Jf j \JY m + Xm * f)m (not relabeled) and Wf G L2(M 2 ;M 2 ) such that 



10/ in L 2 ( 



r>2. Tn,2\ 



y/Ym + Xm * / 

A simple consequence of (3.38), (3.39), and (3.40) is that 

r K||2<W 2 (/,/o). 



(3.40) 
(3.41) 
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In addition, since (\/Y m + Xm * T)m converges towards \f~f uniformly on compact subsets of R 2 , we 
readily deduce from (3.40) that 



jf = a/J Wf a.e. in R 2 
Next, since / = a.e. in R 2 \ V s , it follows from (3.39) that 

li/l 2 



I jf I dx 



< 



< 



Y m + Xm * f 
2 

Jf 



(Ym + Xm* f ~ f) dx 



s/Ym + Xm* f 

wi(fj ) 



(Y m + ||Xm * / - /Hoc) 



(1+Y m f Y m — > 0, 



whence, additionally to (3.37), 

j/ = a.e. in R 2 \V f . 
Finally, owing to (3.40) and (3.43), we have 



/ \wf\ 2 dx = lim / 



(3.42) 



(3.43) 



w f ■ ^ dx = , 

\VS V Y m + Xm * J 

and thus wj = a.e. in R 2 \ . This completes the proof of Lemma 3.2 for /. The statements 
(3.33) and (3.34) for g are proved by similar arguments. □ 



4. Convergence of the time discretization 

We pick now r > and (fo,9o) G ^2- For each integer n > 1, we define (/™ +1 ,<7™ +1 ) G /C2 as a 
solution to the minimization problem 

inf T?(u,v), 
(u,v)eK.2 

where := (/o,3o) and 

«) := i- (Wf («, / T n ) + B W 2 (v, g?)) + S(u, v) , (u, «) G K 2 . 

Recall that (/" +1 ,< +1 ) is well-defined and belongs to # 2 (R 2 ;R 2 ) for all n > 1 by Lemma 2.2. We 
next let (f T ,g T ) : [0,oo) x R 2 — > AC 2 be the function obtained by the method of piecewise constant 
interpolation in AC 2 as follows: (f T ,gr)(t) '■= f° r au * ^ [ nr i ( n + l) r ) an d ^ G N. 

The next lemma collects estimates which allow us to perform the limit r — > and construct in 
this way a weak solution of (1.1). 
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Lemma 4.1. There is C3 > depending only on A, B, a, b, c, /o, and go such that, for all T > 
and t G (0, 1), we have 



(0 

(it) 

(Hi) 
(iv) 

(v) 
(vi) 



\\f T (T)\\ 1 = \\g T (T)\\ 1 = l, 

00 

Y,[wKfr + \fr) + Wl(g n T + \g n T )} <C 3 T, 

n=0 

S(f T (T),g T (T))<£(f ,go), 
[ [f T (T, x) + g T (T, x)} (1 + \x\ 2 ) dx < C 3 (1 + T), 

[\\Af T (s)\\l + \\Ag T (s)\\ 2 2 ] ds<C 3 (1 + T), 



max {T,t} 



k/rll2+ IKrll2 



where 



and 



W fr 



% (-VA(A/ r + 5 5r ) + V(a/ T + bg T )) a.e. in 



a.e. in 



\ V fr 



fg- T (-VA(/ T + 5 r) +cV(/ r + 5r )) a.e. in 7^ 



a.e. in M 2 \ V 9t 



(4.1) 

(4.2) 

(4.3) 
(4.4) 

(4.5) 
(4.6) 

(4.7) 
(4.8) 



Proof. The assertion (4.1) follows from the fact that £ /C2 for all n G N and r > 0. We next 

observe that, since g™) > (/™ +1 , for all n G N, we have 



1 

27 



[Wi(f? + \f?) + B W%(f? +1 ,g?j\ + £(f?+\g? +1 )<£(f?,g?) 



and therefore, for all N G N, 
AT-l 

2r 



1 w- 1 

7- e [^(rvn+B^r 1 .^)] +eu?^)<£u^)- 



(4.9) 



n=0 



Recalling that the functional £ is bounded from below by Lemma 2.1, we obtain (4.2) after letting 
N ->■ 00 in (4.9). Moreover, given T > 0, we choose N > 1 such that T G [JVr, (TV + l)r) in (4.9) 
and arrive at (4.3). Next, the bound (4.4) follows readily from (4.2) and the property (/o,5o) £ K-2 
in a similar manner as (2.9). 

In order to deduce (4.5), we infer from Lemma 2.2 that, for n G N, 

(A - B)\\Af?+ l \\l + B\\A(f T n+1 + g n T +1 )\\l + (a - 6) + ||V/ r n+1 ||l + 6||V(/ T " +1 + g n T +1 )\\ 2 2 
< - [H(f?) - H{f? +1 ) + B (H{g n r ) - H{g n T +1 ))} + (b - a) + ||V/- +1 ||i 
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Summation from n = to n = N— 1 yields 

(N+1)t 

[(A-B)\\Af T (s)\\ 2 + B\\A(f T + g T )(s)\\ 2 

+ (a - 6)+||V/ T (*)||! + b ||V(/ T + g T )(s)\\l) ds 
< [H(f ) - H(f T N ) + B (H(g ) - H(g?))] + (b - a). 



(4.10) 



(N+1)t 



\Nf T (s)\\ 2 2 ds. 



We now use Lemma 2.1, Lemma A. 3, and the estimates (4.3) and (4.4) to obtain 

-(JV+l)r 



[(A - J3)||A/ T ( a )||| + B\\A(f T + g T )(s)\\ 2 2 ] ds 

<C H + [ fo(x) (l + \x\ 2 )dx+\\f \\ 2 2 + C H + [ f?{x) (l + \x\ 2 )dx 
Jr 2 Jr 2 

+ b(c h + [ go{x) (l + \x\ 2 )dx+\\go\\ 2 2 + C H + [ £ (x) (1 + \x\ 2 ) dx 
\ Jr 2 Jr 2 



+ 



4(6 -a). 



(N+1)t 



(f(/r(s),ft-(s)) + Ci) ds<C (1 + T), 



A — B 

for T 6 [JVr, (JV + l)r), hence (4.5). Finally, (4.6) follows from (3.34) and (4.2). 



□ 



Using uniform estimates from Lemma 4.1, we now establish the time equicontinuity of the family 
(f T ,g T ) T . This step is one of the arguments needed to prove the compactness of (f T ,g T ) T . 

Lemma 4.2. There exists a positive constant C4 such that, for all t £ [0, 00), s £ [0, 00), and 
t G (0, 1) we have 



\\fr(t) ~ fr(s)\\ H -4 + \\g T (t) - g T (s)\\ H -4 < C^\t-s\+r. 



(4.11) 



Proof. Let < s < t with s G [vt, {v + l)r), v > 0, and i G [Nt, (N + l)r), iV > u, be given. By 

ioo/Tn>2\ 



virtue of (3.1), (3.30), and (3.32) we have for £ G Cg°(R 2 ) and n > 1, 



JR 2 



< T 



j f nV£dx 



\\D\\\ooW 2 {f?,fr- 



+ 



Using (3.34), (4.1), and Holder's inequality, we obtain 



/ (tf-fr^tdx 

Jr 2 



< r\\ml /2 \\wfn\\ 2 nveiioo + ii^eiu wfa^/r 1 ; 

whence, owing to the continuous embedding of H 4 (M. 2 ) in W^R 2 ), 

/ (ff-f?- 1 )^ <c[w 2 {f^fr l )+w 2 {f^fr 1 )] \\s\\h*- 

Jr. 2 
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Ur(t)-fAs))Cdx 



< E / (ft-fr 1 )^ 

n=u+l jR2 

N 

<cu\\m E [^(/^/r^+^K/^/r 1 )] 

( N \ X / 2 

Vn=i/+1 / 



n=f+l 



<c[V(iV-^)r + ^] U\\m<CVt-s + r U\\m 



which yields (4.11) for f T . A similar computation based on (3.2), (3.31), and (3.33) gives (4.11) for 



□ 



We are now in a position to study the compactness properties of (f T ,g T ) T as r — > 0. 



Lemma 4.3. There exist nonnegative functions f and g in C([0, oo), L2(M?)) and a subsequence 
(,T~k)k>i which converges to zero such that, for all t > 0, 



(f Tk (t),g Tk (t)) -> (f(t),g(t)) in L 2 (1R 2 ;M 2 ), 



L 2 ;M 2 ^ 

rl /Tn>2. Tn>2\ 



(4.12) 
(4.13) 



(/r fe ,3rj -> (f,g) in L 2 (0,t;H l (R 2 ;R 2 )), 
and (f(t),g(t)) G /C 2 . Moreover, we /mwe (f,g) G L 2 (0, t; if 2 (M 2 ; M 2 )) and 

(fr h ,gr k )^(f,g) mL 2 (8,t;H 2 (M 2 ;M. 2 )) (4.14) 

for allt > and 5 G (0,i). 

Proof. On the one hand, we remark that (2.5) together with the estimate (4.3) and Lemma 2.1 
imply that 

(/t)t€(o,i) and (5r)re(o,i) are bounded in Loo(0, oo; H l (R 2 ; M 2 )). (4.15) 
By interpolation, we have the inequality 

IN| 2 <ll<fll<- 4 , heH\R 2 ), 

which gives, together with (4.11) and (4.15), 

||/r(*) " fr(s)\\2 + \\9r(t) ~ 9 r(s)h < C(\t - s\ + t) 1 ^ (4.16) 

for all r G (0, 1), t > 0, and s > 0. On the other hand, for each t > 0, the sequence (/r(i))3r(i))re(o,i) 
lies in a compact subset of L 2 (R 2 ;]R 2 ) by (4.4), (4.15), and Lemma A.l. Owing to these two 
properties, we can invoke [2, Proposition 3.3.1] to conclude that there exists a function (/, g) G 
C([0,oo),L 2 (IR 2 ;]R 2 )) and a subsequence G (0,1), — > 0, such that (4.12) holds true. In 
addition, we deduce from (4.12), (4.15), and the Lebesgue dominated convergence theorem that 

(fr k ,g Tk )^(f,g) in L 2 ((0,T) xM 2 ) for all T>0. (4.17) 
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We improve now this convergence. Given t > 1 and 5 £ (0,1), the estimates (4.4), (4.5), and 
(4.15) ensure that 

f [\\fr(s)\\ 2 H 2 + \\g T (s)\\ 2 H 2] ds + sup {/ (f T + g T )(s,x) (1 + M 2 ) dx) < C (1 + t). (4.18) 
J 8 se(S,t) Urn 2 J 

By Lemma A.l, H 2 (R 2 ) n Li(R 2 , (1 + |rr| 2 ) dx) is compactly embedded in ff 1 (IR 2 ), which in turn is 
continuously embedded in L2(R 2 ), and we infer from [20, Lemma 9] that 

(fr k ,9r k ) -> (/,<?) in La^tj^Ott 2 ,^)), 

which can be improved to (4.13) by using (4.15). Observing next that the right-hand side of 
(4.18) does not depend on 5, we realize that it follows from (4.18) that, after possibly extracting 
a subsequence and using a diagonal process, we may assume that (f,g) e L 2 (0,t;H 2 (R 2 ;R 2 )) and 
that (4.14) holds true. 

It remains to check that (f(t),g(t)) belongs to IC2 for all t > 0. Owing to (4.12) and (4.15), we 
readily obtain that f(t) and g(t) are both nonnegative and in ff 1 (M 2 ). In addition, (4.12) and (4.18) 
imply that (f Tk {t), g Tk (t))k>i converges towards (f(t),g(t)) in Li(R 2 ; R 2 ) from which we deduce that 
II /(*) 111 = llfl'COIIl = 1- Using once more (4.18), this convergence also guarantees that both f(t) and 
g(t) belong to Li(lR 2 , (1 + \x\ 2 )dx). Consequently, (f(t),g(t)) € JC2 for all t > and the proof of 
Lemma 4.3 is complete. □ 

Proof of Theorem 1.1. Let us first check that the functions (/, g) constructed in Lemma 4.3 enjoy 
the regularity (i) and (ii) stated in Theorem 1.1. The boundedness and integrability properties (i) 
follow at once from (4.15) and (4.18) by Lemma 4.3. We next use (4.12) to pass to the limit k — > 00 
in (4.16) and obtain 

\\f(t) ~ f(s)h + \\g(t) - g(s)\\ 2 < Clt-sl 1 ' 10 for all (s,t) € [0,oo) 2 , (4.19) 

which gives the assertion (ii) of Theorem 1.1. 

We now identify the equations solved by (f,g). For that purpose, we use relations (3.1) and (3.2) 
to obtain, for N > 1, t € [Nt, (N + l)r), and f G Qf 3 (M 2 ), 

f f(N+l)r r 

/ (fr(t) -fo)tdx+ / [A(Af T + Bg T ) div(/ T V0 + f T V(a/ T + bg T ) ■ V^] dx ds 

Jm? Jt Jr 2 

linVll N 

< E^^/r 1 ) (4.20) 

1 n=l 

and 

/• /■(A r +l)r r 

/ 07r(t) -g )£dx+ / [A(/ T + (?r) div( 5r V0 + c g T V(/ t + ff r ) • V^] dx ds 

it 2 Jt JR 2 

^^E^Ur 1 ). (4.21) 



n=l 
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Let now t > be fixed. Before passing to the limit r — > in (4.20) and (4.21), let us point out 
that, owing to (4.15) and (4.18), we have for all integers v > 1 and r > with [y + l)r < t + 1, 

Ku+ir i ^ ^ div(/ r V£) + / T V(a/ T + 6 5 r) • V£] 



/ ( ' +1)T [(,4||A/ T || 2 + i?||A ffT || 2 )||/ r ||^ Heilws, 

+ H/rlb HVeiloo (a||V/ T || 2 + 6||V 5r || 2 )] rfs 

<c7(i+t)v^nefe 



< 



(4.22) 



and 



/ ( ] I [A(/ T + 5t ) div( 5r Ve) + c 5t V(/ t + g T ) ■ V£] cfe 

Jut JR 2 



< 



[(||A/ t || 2 + ||A 5t || 2 )|| 5t || h i iieii^ 

+ c||<7r|| 2 IIVCIU (||V/ T || 2 + ||V 5t || 2 )] ds 

<c(i+t) V?U\\w*- 



(4.23) 



We fix 5 £ (0,t). For each k > 1, there are integers and such that t £ [iVfeTfe, (iV/% + 1)t^) 
and 5 G [i/jfcTfc, + l)r fc ). In virtue of (4.2), (4.20), and (4.22) we obtain that, for £ G Q,°(R 2 ), 



+ f I [A(Af Tk +Bg Tk )dw(f Tk VO + fr k V(afr k +bg Tk )-V^]dxds 

oo 



p=l 

-(^ fe +l)r fc 



/ / [A(A/ Tfc + Bg Tk ) div(/ Tfc V£) + / Tfc V(af Tk + bg Tk ) • V£] dxds 

J8 il 2 



+ 



(iV fc +l)r fe 



[A(^/ r , + 5^) div(/ Tfc V£) + f Tk V(a/ Tfc + 6 5rfc ) • V£] dxda 



<c 3 r fc n^eiu + c(i + 1) uiiwi 

<c(i + t) uww* jr k . 



(4.24) 



Let us now pass to the limit — > in (4.24). We note that the convergences (4.13) and (4.14) 
guarantee that 



A(Af Tk + Bg Tk ) Vf Tk - A(Af + Bg) Vf in ^((tf, t) x R , 
A(,4/ Tfc + B ffr , ) /_ -^A(Af + Bg) f in L x ((5, t) x R 2 



2. TID 2 \ 



(4.25) 
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while (4.13) implies 

f Tk V(af Tk + bg Tk ) — > / V(a/ + bg) in Li((0,i) x R 2 ;R 2 ). (4.26) 
We then let k oo in (4.24) and use (4.12), (4.25), and (4.26) to conclude that 

/ (/(*) " /(<*)) e dx + /* / [A(A/ + 5 5 ) div(/V0 + / V(o/ + 65) • V£] dz = 

for all £ G Co°(]R 2 ). By Lebesgue's dominated convergence theorem and (4.19) we may let 6 — > 
and thus obtain the first identity of (1.7). The second identity of (1.7) follows in a similar way, 
starting from (4.21) and (4.23). 

Let us now prove (1.8). We fix t > 0, 5 G (0,t) and take k > 1 sufficiently large so that < 5 
and rifc > 1 such that i G [^fc T fc> («& + l)Tfc). It follows from (4.10) and (4.15) that 

D H (f Tk (s),g Tk (s))ds< / [D H (f Tk (s),g Tk (s)) + (b-a) ||V/ Tfc ( S )|| 2 ] ds 

+ (a-b) f\\Vf Tk (s)\\lds 
<H(f ) - H(f Tk (t)) + B (H(g ) - H{g Tk (t))) ( 4 - 27 ) 



T k (s)\\ 2 2 ds+ / 



+ \b~a\ [ I ||V/ Tfe ( S )|||d S + / \\Vfr k (s)\\ 2 2 ds 



<H(f ) - H(f Tk (t)) + B (ff( 5o ) - # (5r fc (*))) + C \b - a 5 . 
Now, on the one hand, we infer from (4.13) and (4.14) that 

liminf f D H (f Tk (s),g Tk ( S ))ds =limmf f [D H (f Tk (s),g Tk (s)) + (6 - a) ||V/ r j!] ds 

+ hm(a-6) f\\Vf Tk {s)\\ 2 2 ds 

fc->oo J s 

> f [D H (f(s),g(s)) + (b-a) ||V/|||] ds 



+ (a-b) / \\Vf(s)\\ 2 ds = / £>*(/(*), </(*))<*«. 

On the other hand, it follows from (4.4), (4.12), and (4.15) by classical arguments that 

lim H(f Tk (t)) + BH(g Tk (t)) = H(f(t)) + BH(g(t)) . 

k— >oo 

see [13] for instance. Thanks to these two properties, we can pass to the limit k — > 00 in (4.27) and 
obtain 



J D H (f(s),g(s)) ds < H(f ) - H(f(t)) + B (H(g ) - H{g(t))) + C \b - a\ 5 



(4.28) 



for all 5 < t. By the monotone convergence theorem and the assertion (i) of Theorem 1.1 we may 
let 6 ->• in (4.28) and end up with (1.8). 
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In order to obtain the last estimate (1.9), we deduce from (3.34) and (4.9) that, if t > 5 > and 
k is sufficiently large (so that < 5), then 







2 






2S(f ,g )> 


W fr k 


+ B 

2 


W 9r k 


:) 



ds + 2£(f Tk (t),g Tk (t)) 



(4.29) 



the functions Wf T and w gTk being defined in (4.7) and (4.8), respectively. Since £ is bounded from 
below by Lemma 2.1, we infer from (4.29) that {wf Tk )k and (w gTk )k are bounded in L 2 ((<5, oo) x 
M. 2 ; IR 2 ) for all 5 > 0. Therefore, after possibly extracting a subsequence and using a diagonal process, 
we find vector fields Vf and V g in L 2 ((0, 00 ) x M 2 ;M 2 ) such that 



W/rfc ,^ Tfc ) - (Vf , Vg,) in L 2 ((5,oo) x R 2 ;R 2 ) 



for all 5 > 0. 



(4.30) 



Owing to (4.12), (4.15), and (4.30), we can first perform the liminf k — > oo in (4.29), then take the 
limit as 5 — > with the help of the monotone convergence theorem in the resulting inequality, and 
thus arrive at 

2£(f ,go)> [\\\V f f 2 + B\\V g f 2 ) ds + 2£(f(t),g{t)) for all t>0. (4.31) 
J o 

It remains to identify the terms Vf and V g . To this end, we remark first that (4.15) ensures that 
(V 'fi~ k )k an d {\/gr k )k are bounded in Loo(0, oo; L 4 (M 2 ; M 2 )), which implies, together with (4.29), 
that the sequences (j /rfe ) k and (j 3Tfc ) k defined by jf Tk = ^JJ^ k w frk and j 9Tfc = JgT k w gTk , k > 1, are 
bounded in L 2 (5, oo; L 4/3 (IR 2 ; M 2 )) for all 5 > 0. Since L 2 (<5, oo; L 4 / 3 (M 2 )) is a reflexive space, there 
are vector fields If and I g in £2(0, 00; L 4/ / 3 (]R 2 ; M 2 )) and a subsequence of (Tf,)k ( n °t relabeled) such 
that 

(l/,/ g ) in L 2 (<5,oo;L 4/3 (lR 2 ;IR 2 )) for all 5 > 0. (4.32) 



(4.33) 



Combining (4.12), (4.30), and (4.32) gives 

I f = v ffVf and I g = ^fg-Vg a.e. in (0, 00) x M 2 



Consider now a test function S € Cg°((0, 00) x IR 2 ; M 2 ). For each k > 1 and n > 1, we choose as 
test function 

S(s, x)ds , x € M 2 , 

'n-r fc 

in (3.32) for (f™ k ,g™ k ) and find, since (f Tk ,gT k ) is constant on [nrt, (n + l)Tfc): 



J nr k 



(n+l)r k 



(n+l)r k 



[ [A(Af Tk + Bg Tk ) div(/ Tfc ~) + f Tk V(af Tk + 6^) • 5] dxds 
Jw, 2 

/ ;7 -Edxds. 

JR 2 

Summing up the previous identity with respect to n > 1 gives 

/"OO /■ /'OO /■ 

/ / [A(^/ Tfc + 5^) div(/ Tfe H) + f Tk V(af Tk + 6 5t J ■Z]dxds= / j> • Hdxds 

J0 JM 2 JO Jr 2 
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for k large enough (such that supp(H) C (r k ,oo) x R 2 ). Due to (4.13), (4.14), and (4.32), we can 
pass to the limit as k — > oo in the above equality and find 

roo p poo r 

/ / [A(Af + Bg)dw(fZ) + fV(af + bg)-Z]dxds = / I f -Edxds, 
JO JR 2 JO JR 2 

that is, 

I f = -V(fA(Af + Bg)) + A(Af + Bg)Vf + fV(af + bg) in P'((0, oo) x R 2 ;R 2 ). (4.34) 

A similar argument allows us to deduce from (3.33), (4.13), (4.14), and (4.32) that 

I g = -V(gA(f + g)) + A(/ + g)Vg + cgV(f + g) in V'((0, oo) x R 2 ; R 2 ). (4.35) 

Collecting (4.31), (4.33), (4.34), and (4.35) gives the last assertion of Theorem 1.1 and completes 
its proof. 

□ 

Appendix A. Auxiliary results 

It is well-known that ff 1 (M 2 ) is not compactly embedded in L2(R 2 ) due to the non-compactness 
of R 2 but that compactness can be restored by an additional decay at infinity as in the following 
lemma: 

Lemma A.l. The spaces H 1 ^ 2 ) n Li(R 2 , \x\ 2 dx) and H 2 (R 2 ) n Li(M 2 , \x\ 2 dx) are compactly 
embedded in L2(R 2 ) andH l {R?), respectively. 

Proof. Let (h k ) k >i be a bounded sequence in H 1 (R 2 ) n Li(R 2 , |x| 2 dx). Without loss of generality, 
we may assume that there is a function h £ H 1 (R 2 ) such that h k — ^ h in H 1 ^ 2 ). Furthermore, the 
Rellich theorem guarantees that (hk\B N )k>i is relatively compact in ^(Bat) for all integers N > 1, 
where is the open disc centered in zero and of radius N and h^jj, N the restriction of h k to 
Btv- We may then extract a subsequence, denoted again by (h k ) k >i, such that (h k ) k >i converges 
(strongly) towards h in L2(Ojv) for all N > 1. 
First, for each N > 1, we have that 

sup < / \hk(x)\ \x\ 2 dx> > lim / \h k (x)\ \x\ 2 dx = / \h(x)\ \x\ 2 dx, 
fe>l [ii 2 J k ^°°Jo N Jo N 

so that /i G Li(R 2 , \x\ 2 dx). Next, we have 

\\h k - h\\ 2 2 < [ \h k (x) - h(x)\ 2 dx + / \(h k - h)(x)\ 2 dx 
Jd n 



<[ \{h k - h){x)\ 2 dx + -J- / \{h k -h)(x)\ 2 \x\dx 

Jd n iV Jr 2 \o n 

< [ \(h k - h)(x)\ 2 dx 
Jo N 
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Since H 1 ^ 2 ) is continuously embedding in L^(M. 2 ) and the function h belongs to H 1 ^ 2 ) and 
Li(IR 2 , \x\ 2 dx), there exists a constant C such that 

\\hk ~ h\\l < [ \(h k - h)(x)\ 2 dx + — for all k > 1 and N > 1. 

Letting first k — > oo and then N — > oo, we conclude that (hk)k>i converges towards h in L2(M 2 ) 
and thus that H 1 ^ 2 ) n L 1 (M 2 , \x\ 2 dx) is compactly embedded in L 2 (M 2 ). 

Consider now a bounded sequence (hk)k>i in ff 2 (R 2 ) n Li(M 2 , \x\ 2 dx). Owing to the previ- 
ous result, there exist a subsequence, denoted again by (/ifc)fc>i, and a function h £ i/ 2 (M 2 ) H 
Li(IR 2 , |x| 2 dx) such that (hk)k>i converges towards h strongly in L2(M 2 ) and weakly in H 2 (W 2 ). 
Since 

\\w\\ H i < C WwW 1 ^ \\w\\y 2 , w £ # 2 (IR 2 ) , 

a simple interpolation argument then gives that (hk)k>i converges towards h strongly in H 1 ^ 2 ) 
and completes the proof. □ 

The next result was used in the identification of the Euler-Lagrange equation for the minimizers 
of the minimization problem (2.1). 

Lemma A. 2. Consider h £ H 1 ^ 2 ), ( = (C11C2) £ Cg°(M 2 ;M 2 ) and, for e small enough, define 
Ce := id+e£, je := det(-D£ e ), an ^ h £ := {h Ce^/tie C^ 1 )- Then, there is > suc/i i/iat 
((/t e — /i)/e) e6 (o j£c ) «s bounded in L2(M 2 ). 

Proof. Let us first consider the case /i € C^°(1R 2 ). Then, for e small enough, 

\\he-h\\l= \h £ o C £ - ho Q £ \ 2 j £ dx 
Jr 2 

< 2 / |/f (1 ~ je)2 dx + 2 

Jr 2 Je 
Recalling that j £ = 1 + e div£ + o(e), we have 

\<je<2 and < 2 HCII^ e (A.l) 

for e small enough, and we realize that 

2 / \h\ 2 dx < 16 HCII^ \\h\\ 2 e 2 . 

Jr 2 j £ °° 





71 h t 


2 


/ 


— - h Q 


je 




Je 




l/l 


- ho( £ \ 2 j £ dx 



Next, using once more (A.l), 

2 / \h- ho Q £ \ 2 j £ dx <4e 2 I [I \Vh(x + es({x)) ■ ((x)\ ds ) dx 



1 n 2 



'R 2 VO 

<4e 2 [ \((x)\ 2 f \V(ho( S£ )( x )\ 2 dsdx 







<4^ 2 HCIIL / / :'-> dyds 

JO Jr 2 Jes ° Ces 

<8e 2 HWl \\Vh\\ 2 



\Vh\ 



2 
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for sufficiently small e. Combining the above inequalities gives the claimed boundedness of {h e — h)/e 
m L 2 (M. 2 ) for e small enough. 

The general case h G ff 1 (]R 2 ) next follows by a density argument. □ 

We finally recall some well-known estimates for the functional H defined in (1.4), see, e.g., [13, 
Lemma A.l]. 

Lemma A. 3. Let h be a nonnegative function in Li(M 2 , (1 + x 2 ) dx) n L2(M 2 ). Then hhih € Li(IR 2 ) 
and there exists a positive constant Ch such that 

[ h(x) | In h(x) | dx < C H + f h(x){l + \x\ 2 ) dx + \\h\\l, (A.2) 
Jm. 2 Jr 2 

H(h)>-C H - f h(x)(l + \x\ 2 )dx. (A.3) 
Jr 2 
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